The median recoil direction as a WIMP directional detection signal 
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Direct detection experiments have reached the sensitivity to detect dark matter WIMPs. Demon- 
strating that a putative signal is due to WIMPs, and not backgrounds, is a major challenge however. 
The direction dependence of the WIMP scattering rate provides a potential WIMP 'smoking gun'. 
If the WIMP distribution is predominantly smooth, the Galactic recoil distribution is peaked in the 
direction opposite to the direction of Solar motion. Previous studies have found that, for an ideal 
detector, of order 10 WIMP events would be sufficient to reject isotropy, and rule out an isotropic 
background. We examine how the median recoil direction could be used to confirm the WIMP origin 
of an anisotropic recoil signal. Specifically we determine the number of events required to confirm 
the direction of solar motion as the median inverse recoil direction at 95% confidence. We find that 
for zero background 31 events are required, a factor of ~ 2 more than are required to simply reject 
isotropy. We also investigate the effect of a non-zero isotropic background. As the background rate 
is increased the number of events required increases, initially fairly gradually and then more rapidly, 
once the signal becomes subdominant. We also discuss the effect of features in the speed distribution 
at large speeds, as found in recent high resolution simulations, on the median recoil direction. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs), and 
in particular the lightest neutralino in supersymmetric 
models, are a well motivated dark matter candidate [J- 
3j. WIMPs can be detected directly, in the lab, via the 
elastic scattering of WIMPs on detector nuclei [ij . Ex- 
periments now have the sensitivity required to probe the 
theoretically favoured regions of parameter space 043 
and the CDMS experiment has recently observed two 
events in its WIMP signal region [5( . 

Neutrons, from cosmic-ray induced muons or natural 
radioactivity, can produce nuclear recoils which (on an 
event by event basis) are indistinguishable from WIMP 
induced recoils. Furthermore perfect rejection of other 
backgrounds is impossible, for instance 'surface events' 
(electron recoils close to the detector surface) in the 
case of CDMS. As highlighted by the CDMS events, as 
well as the long-standing DAMA annual modulation sig- 
nal Q , demonstrating the WIMP origin of a putative sig- 
nal is absolutely crucial. The direction dependence of the 
WIMP scattering rate (due to the Earth's motion with 
respect to the Galactic rest frame) [9j provides a potential 
WIMP 'smoking gun'. Assuming that the WIMP distri- 
bution is predominantly smooth, the peak WIMP flux 
comes from the direction of solar motion (towards the 
constellation Cygnus) and the recoil rate is then peaked 
in the opposite direction. The recoil rate, in the Galactic 
rest frame, is highly anisotropic; the rate in the forward 
direction is roughly an order of magnitude larger than 
that in the backward direction. A detector capable of 
measuring the nuclear recoil vectors (including the sense 
+p versus -p) in 3-dimensions, with good angular res- 
olution, could reject isotropy of the recoils with only of 
order 10 events [HHHI. Most, but not all, backgrounds 
would produce an isotropic Galactic recoil distribution 



(due to the complicated motion of the Earth with respect 
to the Galactic rest frame). An anisotropic Galactic re- 
coil distribution would therefore provide strong, but not 
conclusive, evidence for a Galactic origin of the recoils. 



Confirmation of the WIMP origin could be obtained by 
verifying that the properties of the anisotropy match the 
expectations for WIMP induced recoils. Assuming the 
WIMP distribution is dominated by a smooth compo- 
nent, the median inverse recoil direction should be com- 
patible with the direction of solar motion. 



To summarise, a WIMP search strategy with a direc- 
tional detector could be divided into a sequence of phases: 

1. Search phase (detection of non-zero recoil signal) 

2. Detection of anisotropy 

3. Study of properties of anisotropy 

which require successively larger numbers of events (and 
hence larger exposures). 



The initial simple search phase aims to detect an 
anomalous recoil signal above that expected from back- 
grounds. To claim an anomalous signal inconsistent with 
zero at 95% confidence requires 4-5 events. The second 
step, as discussed above, would be to check whether the 
Galactic recoil directions are anisotropic and would re- 
quire of order 10 events (assuming zero background). 
In this paper we focus on the third phase, examin- 
ing how measuring the median recoil direction could be 
used to provide confirmation of the WIMP origin of an 
anisotropic recoil signal. In section |TT] we describe the 
input to our Monte Carlo simulations. In section [TTT1 we 
present our results, before concluding with discussion in 
Sec EH 
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II. MODELLING 

We use the same statistical techniques and methods 
for calculating the directional nuclear recoil spectrum as 
in Refs. [llT[l3l ]. We briefly summarise these procedures 
here, for further details see these references and Ref. (l4| . 

A. Detector 

Many of the directional detectors currently under de- 
velopment [HI HH are low pressure gas time projec tion 
chambers (TPCs), e.g. DMTPC [17], DRIFT [SGil and 
NEWAGE 20]. We therefore simulate a fairly generic 
TPC based detector. We assume that the the recoil di- 
rections, including their senses, are reconstructed per- 
fectly in 3d. These are optimistic assumptions, there- 
fore our results provide a lower limit on the number of 
events/exposure required by a real TPC based detector. 
For concreteness we use a S target with an energy thresh- 
old of 20 keV. Finite angular resolution does not signif- 
icantly affect the number of events required to detect 
anisotropy [ill . l2ll ] , provided it is not worse than of order 
10°. Axial and/or 2-d read-out would, however, signifi- 
cantly degrade the detector capability [Ill4l3l |21| . 



the deviation of the peak WIMP flux. It will also de- 
pend on the (a priori unknown) velocity and density of 
the feature. For instance, for a stream with velocity, in 
Galactic coordinates, v str = (—65, 135, — 249) km s -1 the 
difference between the inverse median direction and the 
solar direction only exceeds 5° if the fraction of the local 
density in the stream exceeds 10% [ll| . 

Since the effect of a feature in the speed distribution 
on the median recoil direction is expected to be small 
we do not investigate it in this study. If, with a large 
number of events, a statistically significant deviation of 
the inverse median direction were found, its origin could 
be investigated by studying the energy dependence of 
the deviation. We defer an investigation of this to future 
work. 

Finally, direct detection experiments probe the ultra- 
local dark matter distribution on scales many orders of 
magnitude smaller than those resolved by simulations. It 
is not (and may never be) possible to directly calculate, 
or otherwise measure, the dark matter distribution on the 
relevant scales. In this case if WIMPs are detected, then 
the directional recoil rate could be used to reconstruct 
the ultra-local dark matter distribution [251. 



III. RESULTS 



B. WIMP properties and distribution 

The detailed angular dependence of the recoil rate 
depends on the exact form of the WIMP velocity dis- 
tribution pU [HI HI- However, if the WIMP veloc- 
ity distribution is dominated by a smooth component 
the main features of the recoil distribution (rear-front 
asymmetry, median direction opposite to the direction 
of solar motion) are robust (see e.g. Ref. [ll|)). For 
concreteness we use the standard halo model halo, an 
isotropic sphere with local density p — 0.3 GcVcm~ 3 and 
a Maxwellian/gaussian velocity distribution with three 
dimensional velocity dispersion a v = 270kms _1 , and fix 
the WIMP mass at m x = 100 GeV. 

Numerical simulations find velocity distributions with 
stochastic features at large speeds [23|, [HI. Kuhlen et 
al 



The recoil rate is peaked in the direction opposite to 
the direction of solar motion. To allow ease of compari- 
son with the direction of solar motion we use the inverse 
recoil directions (i.e. the directions from which the recoils 
originate) in our analyses. 

We first examine how the accuracy with which the me- 
dian Galactic recoil direction is determined depends on 
the number of WIMP events. The median direction is 
defined as the direction x mc d which minimises the sum 
of the arclengths between x me d and the individual inverse 
recoil directions x; [26] . It is found by minimising 
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(1) 



where N is the number of events. 



13 find that for high speed WIMPs (v > SOOkms^ 1 ) For a given number of WIMP events, iV wimp , we simu- 



the direction in which the WIMP flux is largest can de- 
viate by more than 10° from the direction of solar mo- 
tion. The effect on the median recoil direction will be 
substantially smaller than this however. Firstly the set- 
up we are considering (S target with energy threshold of 
20 keV and m x = 100 GeV) is sensitive to WIMPs with 
v > 225 kms -1 . Therefore a high speed feature will only 
contribute a small fraction of the WIMP flux. This will 
be true in general unless the WIMP mass is small and / or 
the threshold energy is large. Secondly, see e.g. fig. 3 
of Ref. [IlJ, due to the elastic scattering the recoil rate 
is less anisotropic than the WIMP flux. The deviation 
of the inverse median direction from the direction of so- 
lar motion will therefore be substantially smaller than 



late 10 4 experiments and in each determine the direction, 
Xmcd, by minimizing Eq. ([TJ and hence calculate A, 



A = cos 1 (x me d-x ) , 



(2) 



the angle between the median direction and the direc- 
tion of solar motion, x Q . In Fig. 1 we plot the 50% 
and 95% percentiles of the distribution of A as a func- 
tion of N, 



wimp 



We also investigate the effects of non-zero 
(isotropic) backgrounds. We parametrize the background 
rate in terms of the fraction of events which are signal, 
A = S/(S + B) where S and B are the signal and back- 
ground rates respectively (c.f. Ref. [27]). For compari- 
son we also plot the 5% and 50% percentiles for an purely 
isotropic distribution. 
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FIG. 1: The distribution of the angle between the me- 
dian inverse recoil direction and the direction of solar mo- 
tion, A, as a function of the number of WIMP events, 
A'wimp, for varying signal fractions (from top to bottom) 
A = 0.0625, 0.125, 0.25, 0.5 and 1. The solid and dashed lines 
are the 95% and 50% percentiles. The dotted lines are (from 
top to bottom) the 50% and 5% percentiles for an isotropic 
recoil distribution. 



We now determine the number of events required to 
confirm the direction of solar motion as the median in- 
verse recoil direction at 95% confidence. We do this using 
the same methodology as in Ref. [ll[ ■ Briefly, we use the 
distribution of A for WIMP induced recoils and for the 
null hypothesis of isotropic recoils, to calculate the re- 
jection and acceptance factors, R and A. The rejection 
factor gives the confidence level with which the null hy- 
pothesis can be rejected given a particular value of the 
test statistic, while the acceptance factor is the probabil- 
ity of measuring a larger value of the test statistic if the 
alternative hypothesis is true, We then find the number 
of WIMP events required to give A = R = 0.95 i.e. to 
reject the median direction being random at 95% confi- 
dence in 95% of experiments. This is shown in Fig. [21 as 
a function of the signal fraction A. For zero background 
(i.e. A = 1) 31 events are required, a factor of ~ 2 more 
than are required to simply reject isotropy. The number 
of events required increases as A is decreased, initially 
fairly gradually and then, once the signal becomes sub- 
dominant (A < 0.5), more rapidly. 

Billard et al. [27] have recently investigated using a 
sky-map based likelihood analysis to probe the correla- 
tion of the directional recoil rate with the direction of 
solar motion (and hence confirm the WIMP origin of a 
signal). Their results are qualitatively consistent with 
ours. For instance they find that, for a signal fraction 
A = 0.5, the peak signal direction can be confirmed to 
be within 20° of the direction of solar motion with as 
few as 25 WIMP events. The median direction method 
is however faster and more robust. It particular it has 
the advantage of being model independent (i.e. one does 



FIG. 2: The dependence of the number of WIMP events, 
A'wimp, required to reject the median direction being random 
at 95% confidence in 95% of experiments on the signal frac- 
tion, A. 



not need to assume an exact form of the WIMP velocity 
distribution) . 

If desired, further confirmation of the WIMP origin of 
an anisotropic recoil signal could be obtained by studying 
the variation of the median recoil direction in the lab 
frame 18]. Over the course of a sidereal day the peak in 
the recoil distribution traces out a small circle on the sky. 
For a detector located in the Northern hemisphere during 
a sidereal day the peak recoil direction in the lab rotates 
(roughly) from down to South and back. A periodogram 
analysis could be used to verify that median direction 
varies over a sidereal, rather than solar, day (c.f. Ref. [28[ 
for the annual modulation). 



IV. DISCUSSION 

In this paper we have investigated how the median re- 
coil direction in directional detection experiments can be 
used as a WIMP signal. Assuming a smooth WIMP dis- 
tribution, the peak WIMP flux is from the direction of 
solar motion and the median recoil direction is in the op- 
posite direction. An ideal detector could reject isotropy 
of recoils with only of order 10 events [13, [ll|. Con- 
firmation of the WIMP origin of an anisotropic recoil 
signal could be obtained by studying the details of the 
anisotropy and in particularly confirming that the me- 
dian inverse Galactic recoil direction coincides with the 
direction of Solar motion. We find that, with an ideal 
detector and zero-background, to confirm the direction 
of solar motion as the median inverse recoil direction at 
95% confidence requires 31 events (see also Ref. j27|). 
Non-zero isotropic background would increase this num- 
ber, significantly if the signal is sub-dominant. 

For concreteness (and in the absence of a definitive al- 
ternative) we have used the standard halo model halo, 
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which has an isotropic Maxwellian speed distribution. 
While the detailed angular dependence of the recoil rate 
depends on the exact form of the WIMP velocity distri- 
bution, if the WIMP velocity distribution is dominated 
by a smooth component the median inverse recoil dis- 
tribution will be close to the direction of solar motion 
(see e.g. Ref. [ll|). Recent high resolution simulations 
have found stochastic, features in the speed distribution 
at large speeds [H, |2J] . The direction in which the flux 
of high speed WIMPs is largest can deviate by more than 
10° from the direction of solar motion 24J. However the 
deviation of the median inverse direction will be small 
compared with that expected from an isotropic recoil di- 



rection. We therefore conclude that features in the speed 
distribution at high speed are unlikely to affect the util- 
ity of the median recoil direction as a WIMP signal. If, 
with a large number of events, a statistically significant 
deviation of the inverse median direction were found, its 
origin could then be investigated by studying the energy 
dependence of the deviation. 
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